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^ Abstract 

' In the mean-field regime, we prove convergence (with exphcit bounds) of the 

many-body von Neumann dynamics with bounded interactions to the Hartree- 
. von Neumann dynamics. 



1 Introduction 



Derivation of macroscopic equations from microscopic ones is one of the main 
I challenges of Mathematical Physics. This is usually a daunting task met with 

a very limited success. In the last few years a considerable progress was made 
on one such problem - derivation of the Hartree, Hartree-Fock and Gross- 
Pitaevskii equations in the mean-field and Gross-Pitaevskii regimes, respec- 



> 

in 

O; tivelv ( [BEGMYl IBGGMl iFGSl iFKSl iFKPl M iLSYl iLSl lErYl lESYl lESl iKFl 

On . IKMIIGM , KSSj). Though the work on the Gross-Pitaevskii limit is quite recent, 

' the work on the mean- field one goes back to the papers \}le\ \GV\ [S] . 

In this note we prove the convergence of solutions of the A^— body von 
^ . Neumann equation with product initial conditions to the A^— fold product of 

solutions of the Hartree-von Neumann equation with the corresponding initial 
conditions and estimate the rate of this convergence. The regime we consider 
is the mean-field one, i.e. with the number of particles going to infinity, while 
the strength of interaction decreasing in the inverse proportion to the number 
of particles. In our analysis we follow closely the beautiful work [FGS| . One of 
the new elements of our approach is a Hamiltonian formulation of the Hartree- 
von Neumann equation. While this work has been written up there appeared 
e-prints [RS| . [ErSj and [GMMj giving, by different techniques, estimates of 
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the rate of convergence in the case of the A^— body Schrodinger and Hartree 
equations. 

Let := p ... (8) p. We start with the time-dependent von Neumann 
equation for a system of N bosons 



at 

PN\t=0 = Po" , 



with pn = /97v(t) acting on L^(]R^)®^ and po, a positive, trace-class operator 
on of trace 1. Here Hn = H% + Vn, with H% = - Xl^i = h is 

a self-adjoint operator in a variable x, e.g. hx ■= + W{x), and 

N 

Since v{xi — Xj) = Yli^j ~ ^j) +v{xj —Xi)) we can assume without 

loss of generality that the two-body potential v is even: v{x) = v{—x). We 
consider the mean field regime: N ^ oo and g ^ with gN c. By changing 
V, if necessary, we can assume that 

We will relate the von Neumann equation ([1]) to the Hartree-von Neumann 
equation 

ih^ = [h + {v*np),p] 
p\t=o = Po, 

where v = v{x) is the same two-body potential as above and np{x, t) := 
p{x;x,t), the probability or charge density, with p{x;y,t) the integral kernel of 
P- 

For V € L°° one can show easily that (jH) is globally well-posed on the space 
of positive, trace-class operators and that the trace, Trp, and the energy, 

E{p) := Tr{hp) + j Up v * Up, (5) 

are conserved. Moreover, p is non- negative, provided so is po- See Appendix 

m 

In Section [4] we show that ([4]) is a Hamiltonian system with the Hamiltonian 
([5]) and the Poisson bracket 

{Aip), B{p)} = -'-Tr {dpA{p)pdpB{p) - dpB{p)pdpA{p)) , (6) 

where A[p) and B[p) are differentiable functionals of p and the operator (Frechet 
derivative) dpA{p) is defined by the equation Tr{dpA{p)^) = dsA{p-\- s^)|s=o- 
This, as was mentioned above, plays an important role in our analysis. 
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Finaly, note that since the integral kernel of the operator /9^|j=q is symmet- 
ric with respect to permutations of particle coordinates, the same is true for 
the solution to ([1]). 

To formulate the main result we need some notation and definitions. For 
any Banach space X, we denote the space of bounded linear operators from 
X to itself by B{X). Let L^R^^^) be the subspace of L'^{R^^^) consisting 
of the functions that are symmetric with respect to permutation of particle 
coordinates, and let 

P|^^'(xi,...,xm) := ^ ^'(x,(i),...,x,(M)), (7) 

where denotes the permutation group of the set {1,2,...,M}, be the or- 
thogonal projection onto the subspace L|(IR^*^) of L^(]R^*^) We denote by I'^ 
the identity operator acting on q coordinates. 

Definition 1. Let Ap := i?(L|(]R'^P)). For p < N, we define the maps (j)^ : 
Ap An by 

<(«):= P#(a®/^-P)P#. (8) 

Let An,p be the image of Ap under (pp . Its elements will be called quan- 
tum p-particle observables or simply p-particle observables. Note that An,i C 
An,2 c ... An,n = An- 

The following theorem relates the asymptotic behavior of ([1]) as — > cxo 
with (H]). 

Theorem 1.1. Assume that v is bounded and even and that 1^ holds. Let pN 
solve ^ with Trpo = 1. Then for any p and any A = (pp (a) G An,p, and for 
N ^ oo, 

Tr{ApN) - Tr{Ap'^^) ^ 0, 
where p solves Moreover, we have the estimate 

\Tr{ApN) - Tr (Ap®^) I < 2^m+^^PN~^^'^ Ma,, (9) 

where r = gji^j;^, [7] is the integer part of ^ and 7(t) := -^r^/ 

Remark 1. The N -bound in this theorem is rather poor, especially compared 
with estimates for the Schrddinger equation mentioned above. It is improved 
somewhat in JAv!^ . where an extension of this result to Coulomb-type potentials 
is also presented. 

From now on we fix the particle number N and sometimes drop the corre- 
sponding index from the notation. For instance, we write Ps for Pg , (f)p for 
(^^, and V for Vn- 

The rest of the paper is devoted to a proof of Theorem lA.li The paper is 
organized as follows. In Section [2] we derive a convenient equation for the map 
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Tt{A) := e ft e n Ae » e ft , which is connected to the l.h.s. of Q. 
To this end we use a decomposition of the commutator with the many-body 
potential into the tree and loop operators, introduced in |FGSj . We use this 
equation in Section [3] in order to approximate Tt by an operator Tj^ whose 
expansion contains only tree operators. In Section H] we discuss the Hamilto- 
nian and Liouvillean formulations of the Hartree-von Neumann equation and 
the Dyson expansion for the latter. We also show that in certain (symbolic) 
representation the tree operators act as Poisson brackets. In Section [5] we prove 
the result of Theorem lA.ll for small times and in Section [6] we use the group 
properties of the von Neumann and Hartree-von Neumann dynamics in order 
to extend the proof to all times. 
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2 Map Ft and its equation 

In this section we derive a convenient equation for the family of operators Tt 
acting on elements of and defined by 

Tt{A) := e~f^ e~f^ Ae~f^ e~f^ . (10) 

This family is related to the l.h.s. of ([9]) as Tr{Ap]^) = Tr{Tt{At)pQ^), where 
Po,Pn are the same as in ([T]) and At denotes the free evolution of A: 

At := e~f^ Ae~n~ . (11) 
Writting r,.(A) as the integral of derivative, we obtain 

Tt{A)=A+'- frriiVr, A])dr. (12) 

A simple analysis shows that the solution of this equation is unique in 
C{[0,cx^),B{An)). 

Now we decompose the commutator on the r.h.s. of (|12p in a convenient 
way (cf. |FGSj ). For any i,j < N,i ^ j denote by V^^ the multiplication 
operator by v{xi — Xj). 
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Proposition 2.1. We have for a £ Ap, 



'-[Vr, (t)p{a)] = Tr{(t>p{a)) + Lr{4>p{a)), (13) 



where Tr-{(j)p{a)) = for p > N and otherwise 



Tr^M) = ^i—P<p^^^{Xp,.{a)), (14) 



with Xp^r '■ Ap — > Ap+i, defined by 



Xp^a) ■.= pPl+'-[Vr^\a^I^]Pl+\ (15) 



for p < N, and 



LriU^)) = 4>p{Yp,M), wzth YpAa) = ^^j;^iPl^[Vr-\a]Pl). (16) 

Proof. Recall the notation Vt = e'~r^ Ve ~ . Using equations Vr = 'l2i<j ^r'' 
and ([3]), we obtain 

1 ^ 

[Vr,A] = -Y.iV;^^A]. (17) 

i<j 

Since the operator Vr = jj' Yli<j is permutationally symmetric, we have 
from dHD that 

'-[Vr,A]='-Ps[Vr,a®l''-P]Ps 

l,N 

= ^PsJ2[V;',a ® /^-^])P5 = Tr{A) + Lr{A), (18) 

i<j 



where 



and 



p N 
i=l j=p+l 



1 ■ ^'^ 



By symmetry we have 



TM) = ^-^^^j^Ps{\\Vr^\a®I^-^])Ps. 
Now, since P^^^Pg = PgP^^^ = Psi we have furthermore 

TriA) = ^PsipPl+'C-ivr-"', a ® I])PI+' ® I^-P-')Ps 
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which gives (|14|) - (|15|) . Similarly, we find 



This, due to (|8|), gives □ 

Remark 2. In general, Tr{4'p{a)) ^ Tr{(j)q{b)), even if cj)p{a) = 4>q{h). Thus, 
e.g. the expression T^A should he understood as Tr(j)p{a) with A = (j)p{a). How- 
ever, our abuse of notation will not cause a confusion. 

Tr and Lr will be called the tree and loop operators, respectively (see |FGSj ). 
Observe that the equation (jl4p implies that 



A G An,p =^ TriA) G An,p+i (20) 
and equation (fT6|) implies that 

A e An,p =^ Lr{A) G An,p. (21) 



Combining equations (I12p and (I13p from above we obtain the following 
equation for Tt 

rt{A)=A + [ rr{Tr{A))dr+ [ Tr{Lr{A)))dr. (22) 
Jo Jo 

Introducing the notation T = F. = {Tt,t > 0) we rewrite this equation in a 
more compact way 

r = I + KT + R{T) (23) 
where ^ 

{KG)t := [ GsTsds (24) 
Jo 

and ^ 

R{T)t= [ T.Lsds. (25) 
Jo 

Equation (|24p defines the operator K : G ^ KG acting on the families G = 
{Gt G B{AN),t> 0}. Clearly, is a bounded operator on C([0,r],5(^Ar)),VT 
0. 

Proposition 2.2. Let K he the operator defined in equation (|24|) . Then I — K 
is invertihle and 

N-p 

yA G An,p, {I - K)-^GA = K^'GA. (26) 
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Proof. Let A E An,p- The definition of the operator K imphes that 

{K^G)tA= I rtGt„Tt,,...Tt,A, (27) 

where J^t = Jq dti J^^ dt2--- /q""^ dtn- (Here Afj is the n-symplex < t„, < 
tn-l < ••■"< ti < t.) 

Equations ([20]) and Tr{A) = VA € An,n imply 

r,„T,„_,...r,,yl = 0, VA€^7V,p, n>N-p. (28) 

Hence, by ([28]), {K''G)tA = ior n > N - p. This gives 

oo N—p 

K'^GA = K'^GA. (29) 

n=l n=l 

On the other hand, {I-K) X^^^i K'^GA = GA and ^^^^ K''{I-K)GA = GA, 
which completes the proof. □ 



3 Approximation of F 

Let r(-f^) := J2n=oK''I, which, due to ([29]), is a finite series on An,p- Equiva- 
lently, we write 



rS^):=^/ d-tTt„...Tt,. (30) 
Proposition 3.1. For t <t := — , we /ia?;e 



||(r,-rf)<^,IU^^^, <2P-2^^. (31) 

Proof. Using Proposition 12.21 and equation (j23p we obtain that 



r-r^ = (/-/f)-ii?(r) = ^/f"i?(r). (32) 

n=0 

Using equations (|25p and (|27p . we find that 

{K-R{T))t = [ d^+H Tt„^,Lt„^,Tt„...Tt,. 

n + l 

Using equation (|14p and (|16p . we obtain that 

Lt„+iTt^-Tt,(l)p{a) 

= 4>p+n{^j^ ^ ^ _-^^^'|^^yp+„,t„+iXp+„_i,f„Xp+„_2,t„_i...Xp,ti (a)). (33) 
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Using equations (fTSj) . ||-P|^|| = 1 and || V^*-' ||b(l2(]r6)) = ||f||oo) we derive the 
estimate on the tree and loop operators 

||^P,t(a)IUp+i < ^^^y^p||«IUp, (34) 
uv ( Ml ^ lbllooP(p- 1) II II ,„^^ 
Equations ([5^ . ([55]) . ([3]) and ||</'p(«)|UiVp ^ II^IUp (since ||-Ps'|| = 1) imply that 

1 1 -C-in+l ^tn • • -^t 1 I Up ^^JV,p + „ 

^ / 2||^|[oo \" (iV-p)! (p + n)!(p + n-l) 
"V ^ y (A^-p-n)!iV" (p-1)! 2Af ^ ^ 

which together with the fact that ||rt||_4^p_^_4^ = 1 and the equality 

(Tt = I dti dt2... / dtn = -, (37) 



A*, Jo Jo Jo 

implies that 

< 1 (±Y''' {N-p)l {p + ny. p + n-1 

-2V4ry {N - p - n)lN'^p - l)l{n + ly. N ' ^ ' 

Furthermore, using the inequalities 

{p + ny. 



(p- l)!(n + l)! 

and 



< (39) 



(N-py , , 

< 1, (40) 



{N -p- ny.N'^ 

we simplify ([3]) as 

n+l 

2r J N 

To conclude our calculations we use the following equality: 



K-R{r))tMA.^^^A. < 2^-^ ( ^i^. (41) 



J]]2""(p + n- 1) = 2p. (42) 



n=0 



To derive (gS}, let A = E^=o 2~"(p+n- 1). Then 4 = E^=o 2"""Hp+?^-1) = 
2~"'(p + n — 2), and therefore subtracting the latter equation from the 
former, we obtain that j = {p - 1) + J2n=i 2"" = P, which gives immediately 
equation (|^ . Equation (|12]) yields for t < r that 

N-p-l 

J2 \\iK^RiT))tMA^,,-.A^<2^"'^-. (43) 

n=0 

This together with (j32|) and ([29j) proves equation (j3T]) for t < t. □ 
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4 Classical Field Theory for the Hartree- 
von Neumann Equation 

In this section we develop Hamiltonian and Liouvillian representations for the 
Hartree equation viewed as a classical field theory. We begin with key defini- 
tions. Let Ii denote the space of all positive, trace class operators on L^(M^). 
For any operator a £ Ap we define the p-particle classical field observable 
^4*^ : Ii — 5- C by the equation 

A^p) := Tr{ap^P). (44) 

We equip the space of these functionals with the norm 

WA^W := supp(zi,^Trp=i\A''{p)\. (45) 

The last two equations imply that 

Ul < \\a\U,. (46) 

For a functional A^(p) we define the operator (Frechet derivative) dpA^{p) 
by the equation Tr{dpA'^{p)^^) = dsA'^{p+ s£,)\s=q. On the space of classical field 
observables we define the Poisson bracket by 

{A'{p\ B'^ip)} = -'-Tr {dpA%p)pdpB-{p) - dpB'{p)pdpA'^{p)) . (47) 

The Jacobi identity is proven in Appendix [Bj Note that if A^ is a p-particle 
classical field observable, and B'^ is a g-particle classical field observable, then 
{A^, B'^} is & p + q — 1-particle classical field observable. 
Furthermore, we observe that 

{A-{p),B%p)}\p=p^ = '-j {d^^,)A-d^^^^B--d^^^^A-d^(^,)B-){^,^)dx, (48) 

where P^, is the rank-one projection on the vector -0 and ^'^('0, V') ■= A'^{P^). 
The r.h.s. is the standard Poisson bracket for the Hartree equation (see |FGS| ). 
Indeed, dpA'^{p) and dpB'^{p) are operators on L^(]R^) (1— particle observables) 
and therefore 

Tr{{dpB''){P^)P^{dpA'){P^)) = {{dpA'){P^ri;,{dpB'){P^)i;). 

Since, as it is easy to see, {dpB'^){P^)il^ = "^^(^.^-^^(V', V') and {dpA'^){Pjf,)*ip = 
a^(,)A^(^,V^), this gives Tr{{dpB^){P^)P^{dpA''){P^)) = / 
which implies the desired relation. 

Introduce the classical Hamiltonian functional 



H^{p) := H''^{p) + V^{p), 



HartreeApprox, April 20, 2009 



where (extending (j44p to the unbounded 1-particle observable h) 

H'^'^ip) := Tr{hp), and V'{p) := ^rr(t>p®2). (49) 

Note that the Hartree-von Neumann equation (j4|) is equivalent to the equation 

dtp = {H%p),p}, (50) 

which motivates the above definition of the Poisson bracket. 

Remark 3. The last equation holds in the weak sense that for all a Ai, 

dtTr{ap) = Tr{a {H^, p}). (51) 

Using the linearity of the Poisson brackets in the second factor, we obtain 

Tr{a{H^,p}) = {H^,Tr{ap)}. (52) 

Thus equation ()50p is equivalent to the equation 

dtTr{ap) = {H',Tr{ap)}, (53) 

or dta^{p) = {H^[p),a^{p)} for all 1 particle observables a. 

To show ()50p we use the definition of the Poisson bracket and the relation 
Tr{dpp£,) = S,, which follows from the definition of dpA'^{p) above to obtain 

{H\ p} = -^ {dpH%p)p - pdpH^ip)) . (54) 
Next, computing dpTT^^p), we conclude that is equal to the r.h.s. of 

(SD- 

Let $j be the flow given by the Hartree-von Neumann initial value problem 
dl]), i.e. $t(/9o) := pt where pt is the solution of (|4|) at time t. We denote by 
the flow of (j3]) for v = (the free flow). We define the Hartree-von Neumann 
evolution on the space of p-particle classical field observables by 

UtiA^p)) := A'^iMp)), (55) 

and the free Hartree-von Neumann evolution by U^{A'^{p)) := A'^{^^(p)). In a 
standard way we derive the following equation 

dtUM'ip)) = Ut{{H'',A'^}{p)) (56) 

(and similarly for Uf{A'^{p))). 

Let Vf^ denote the free evolution, Uf{V'^), of the 2-particle classical observ- 
able V^. A simple computation gives that 

Vfip) := iTriutp""'), (57) 

where Vt is the operator 'ip{xi,X2) eS''^''^i"'"''^2)*7;(xi— X2)e~'s^*^'*"=i+'*"=2)*^(xi, X2). 
In what follows we denote the action of Poisson bracket as 

PvA^ := {V,A''}. (58) 
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Proposition 4.1. We have the following expansion: 

oo 

Ut{A^) = Y,Aln, (59) 



n=0 

where A'j: q = A^ := U^{A'^) and, for n > 1, 

Al^ = [ d^tPv^c Py.A'l, 
with the following the estimates 



(60) 



IKJI<(^)''2P-i||a|U^ (61) 

(in particular, for t < t the series converges in the norm ()45p ). 

Proof. Define A^ := U^^{Ut{A'^)) . By a standard argument we have that dtA^ = 
Pyc^A^. Integrating this equation, we obtain immediately that 



A1 = A'+ dt^Py^^Al. (62) 



Iterating this equation and applying Ui to the result we obtain ([59|) . with 
Al^ := U?{A<^) and := dt, J^*^ dt2 f^^' dt^U^iPyc^^ -Pvf.^^'), n > I, 
which after a change of variables of integration gives (j60p . 

It remains to prove (j6ip . which shows that the series (j59p converges and 
which we are going to prove next. Recall the notation = p ^ ... ® p, the 
A^— fold tensor product. 

Lemma 4.2. Let X^n^f be the operator defined in equation (|15p . Then 

Tr{Xp,{a)p^P-^') = Pv^caip). (63) 

Proof. The proof follows from the relation dpA'^{p) = pTrp-i{aff~^ ® I) ., where 
Trp-i is the partial trace over the first p—1 coordinates, and a simple compu- 
tation. □ 

Iterating ()63|) we obtain equation 

rr(Xp+„_i,t„...Xp,i,(a)/>^^+") = Pv^^^...Pv^^a{p). (64) 

Next, using equation (j34p we obtain that 



The last two equations together with (j45p . the formula = (j)p{at), where 

at = e A ae and the isometry of the free evolution, at, imply 
IIP P ^cii ^ (p + ra- 1)! ^ 2||i;||oo y 

Equations ([371), ([Ml),®, ([Ml) and the definition r := give ([61]), which 

completes the proof of the proposition. □ 
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5 Hartree von Neumann approximation for 

t < T 

In this section we estimate the difference between the quantum A^-body av- 

erage Tr{A{t)p® where A{t) := e ft Ae ~ , and the classical evolution 
Ut{A'^{p)), where p satisfies Trp = 1. 

Proposition 5.1. For all A G An,p o-nd for all t < t we have that 

|Tr(^(i)p®^) - Ut{A^{p))\ < 2^+'P±^i\\a\U,. (67) 

Proof. The proof of this proposition uses the following auxiliary lemma: 
Lemma 5.2. For any p-particle observable A and any n < N — p we have that 

Proof. Let A = 4)p{a) with a € Ap. Then At = (l)p{at), where, recall, at = 
e n ae ~ . Now, using the facts that 



Ttr,---Tt^<t'p{a) - _"^^^",^„ '?!'p+n(^p+n-l,t„^p+n-2,t„_i---^p,ti(Q)), (69) 

which follows from equation ([H]), that Psp®'^ = P®^ and that Tr{^p{a)p'^'^) = 
Tr{ap®P), we find 

{Tt^...TtMt)np) = (^i^:^J),^„ rr(X,+._i,„...X,,,(a,)p^^+"). 

Now, equation (f68ll follows from the last equation and equations ([Mil . and 

At = (l)p{at). □ 

Now, equations (|68p and (|60p imply that for any n < N — p 

This, together with ([28]), (l30|), (04]) and jM]), yields that 

Tr{rf{At)p^'')-UtiA'^ip)) 



N-p 



n=l ^ ^ ^ r / n=N-p+l 

which together with (j45p . Trp = 1 and (I6ip gives for t < r 



(71) 



\Tr{Tf{At)p^'') - Ut{A%p))\ < ^-S^MMa,, (72) 
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where Sp^jy 
form 



En=7 ( 1 



{N~p-ny.N" I 2 



m + X/ 



N-p 



5, 



'p,N '■- 



(N-p)] 



^ (iV-p-n)!iV"2" 
The fohowing inequahty is proven in Appendix [Cl 



oo 1 
n=N-p+l 2" 



2P- 



N-p 



(iV-p)! 1 ^ 2(p+l) 



2P~^ We trans- 
(73) 

(74) 



Equations (l73|) and dS]) imply that 5p,Ar < 2^2+1^ which together with ([72 



and (j3ip impHes that for t < r 



rr(rj(^j)p®^) - C/t(^^(p))j < [^2^^-^ + 2^ 



p 



N 



N ^ 



(75) 



iH]\rt iH^jt 



Recah the notation A{t) = e n Ae ft . Due to the equations (jlOp and 
P!T|) . we have that A(t) = Tt{At). This, together with ([75]) . imphes Proposition 
EH □ 



6 Hartree approximation for arbitrary t 

Now we prove our main result, Theorem lA.ll In what follows A = 4'p{a) is 

a p-particle observable and at{A) = A{t) = e « Ae ~ . We proceed by 
induction. Equations ([67|) and A{t) = at{A) imply that 

|rr(a.(A)p^^) - UM''{P))\ < ^^^^l^llalU,. (76) 

Let Lk = (j^\y^ so that Lfc = kLk+i- We assume that for any A = (pp{a) € An,p 
and for some k > 1 

|rr(afc,(^)p®^) - C/fc,(A^)| < RpMWWa,, (77) 

where 

= 2^=^ (2S'=i-^'-^ + 2-^^) , (78) 

and prove it for k + 1. For A; = 1, (|77p follows from (j76p . since Li > 3 and 
p>l. 

We begin with some preliminary inequalities. Let A = <j)p{a) and 

An{t):= [ Tt„...Tt,At. (79) 



Since At = (f>p{at), where at = e n. ae n is a p-observable, we have by (I69p 
that An{t) := (pp^nio-nit)) with 

{N-py. f 

an{t) = f^jy _ p _ ^yjs^n J Xp+n-l,t„Xp+n~2,t„-i-^P,tii'^t), (80) 
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h 



which, together with (f34l) . (f37l) . (f39l) and (l40]l and the definition r :— g||^y|| — , 
gives 



|A(*)IUiV,p+„ < \\an{tju^p+n 



A.r) 



+ (iV-p)! / 2||^||oor 
^ TnTTTTTT Z i ll"IUp (81) 



{p-l)\nl {N-p-n)im\ h 

. (iV-p)! /^^V 11.11 ^89^ 

(iV-p-n)!iv47j """-^'^ (^2) 



<2P-^""(^-J ||a|U^. (83) 

Using that ar{A) = rT-(^T-) and using ([30]) with ([3T]) and ([83]) . we obtain for 
L<N-p that 

L-1 

||a.(^) - J] A„(r)|U, ^ 2^ + 2'^) ll""-^''- 

n=0 

Next, we claim that for L < N — p 

\Ur{A^) - Ar^irTl < 2M ^ + 2-M ||a|U,. (85) 

n=0 ^ ^ 

Indeed, by (l59|), ([70]) and (IT9|) we have that for L < - p 

L-~l L— 1 (M— V °° 

^.(^^) - 5: A„(r)^ = 5^(1 - + E (86) 

n.=0 n=0 ^ ■'^ ^' n=L 

where Aq{t) = A^. This and equation (|61|) imply 

L-1 

\UMn -^MrYl < Sp,N,L\\a\\A„ (87) 

n=0 

where 5,,^,^ := En=o fl " nvS^lynv^) 2^"'"" + E:r=L2^~""^ Proceeding 



as in Eqns ([73D and ([H]), we obtain Sp^N,L < ^2^ + 2^"-^. This inequahty 



(Af-p-n)!Af"^ 

together with ([87|) gives ([85]). 

Now we prove ()77p for A; + 1 (assuming it for A;). In what follows we use the 
notation (•) = Tr{-p^^). Let s = kr. We have by (j84p and the linearity and 
unitarity of Os that for L^+i < N — p 

n=0 

Next, using this inequality, using that j4„(r) are {p + ?i)-particle observables 
(which follows from equation ([50}) ) and using ([77]) and ([55]) . we obtain 

^fe+ 1 — 1 

|(a.(a,(^)))-[/.( £ ^„(r)'^)| 

n=0 
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-1 



< J2 ^p+n,fc||an(T)|U,+„+2P(^+2-^^+0ll«IUp 
n=0 

Equations ([85]), dMl) and (a^+s(^)) = (as(a,-(A))) imply 

\{ar+s{A)) - Us{Ur{A'))\ < Tp,N\\a\\A„ (89) 

where 



n=0 



We claim that 



Tp,N < Rp,k+i. (90) 



Indeed, 

E i?p+n,fc2^'-i-'^ = 2('=+^)?'+('=-i)"-i(2S-i^'^'-^^ + 2"^^-) 

n=0 n=0 
n=0 

< 2{fc+i)p+(fc-i)ife+i-if2Er=i'-ir-^i±i±ti + 2--^M 
^ iV ^' 

Since = /cLfc_|_i, we find 

n=0 

<; 2('=+i)p-i(2Srii'-^^--f'fe+i-^ + 

This inequality, the definition of Tp^^r and elementary bounds imply the estimate 
(HDD, provided k>2. 

Since s = A;r and since UsiUriA")) = Us+r{A''), Eqns ([89]) and imply 
equation ()77p with /c — > /c + 1. Thus ()77p is shown by induction. Take k — 1 to 
be the integer part of ^, i.e. k — 1 = [^], and let r' := Then t' < t and we 
proceed as above but with r replaced by r' to prove ([77|) with r replaced by 
r'. Next, using that ^^^^ rL, = ^(j^^ r^;^ < 2eLo and taking Lo = ^ 
in (j77p we arrive at 

|Tr(at(A)/)^^) - ^7j(A"(p))| < 2([^]+i)p (^piV-^/^ + N^^^?^'^ \\a\U,. (91) 
Now, by the definition of aj(A) and pN we have that Tr(Qt(^)p®^) = Tr{Apiq), 



and therefore ()9ip implies ([9]). Theorem lA.il is proven. 
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A Appendix: Hartree-von Neumann equa- 
tion 

In this section we scketch proofs of some of key properties of the Hartree-von 
Neumann equation For works on the related Hartree-Fock equation see 
[BDFl [C| ICG| . Let Ii denote the space of all positive, trace class operators on 

L2(M3). 

Theorem A.l. Assume that v is bounded. Then the Hartree-von Neumann 
equation is globally well-posed on Ii and the trace and the energy are con- 
served. 

Sketch of Proof. We will display the i— dependence as a subindex. Let Ji 
denote the space of all trace class operators on L'^{W^). Using the Duhamel 
formula we rewrite (jH) as the fixed-point problem p = F{p) on C([0,r],Ji). 
Here T will be chosen later and 

F{p)t := at{po) -\-i atsHv * npjps)ds. (92) 
Jo 

, , iht iht 1 II II 1 1 II II 

where dtlJ) = e ^ 76 i . Denote the trace norm by ||- ||i and let \\p\\t ■= 
supo<s<T IIPsIIi be the norm on the space C([0,T],Ji). Recall that ||/|[oo 
denotes the L°°-norm of a function /. Let Vx{y) '■= v{x — y). Using that 
\{v * np){x)\ = \Tr{vxp)\ < ||v||oo||yo||i, we obtain 

Wat-siiv * npjps)\\i < ||w * npJIoollPslli < IbHoo ||p||t- (93) 

This estimate shows that F maps any ball in C([0, T], Ii) of radius R > 2 into 
itself, provided T < l/||t;||oo^^- Similarly, one shows that F is a contraction 
on such a ball, if T < l/2||?;||oo-R- Hence our fixed point equation in any 
Br, R > 2, has a unique solution for T = l/2||7;||oo-R^- This solution solves 
also the original intial value problem (jH). 

Since pt and p^ satisfy the same equation ([4]) with the same initial condi- 
tion pq = Pq, we conclude by uniqueness that pt = pi- Since the trace of a 
commutator vanishes, one has that the trace of pt is independent of t. 

We show that the eigenvalues of pt are independent of t. We denote by Aj 
and the eigenvalues and the corresponding eigenfunctions of pt and compute 

5iAfc = dt{(t>k, Pt(t>k) = {Ot4'k,p(t>k) + {4>k, ^[hp,Pt]4>k) + {(t>k,Ptdtct>k)- (94) 
Since p is self-adjoint, pt(l)k = ^k4>k and {(j)k, [hp^, pt]<l)k) = 0, this gives 

dtXk = >^k{dtct>k,4>k) + >^k{4>k,dt(t)k) = hdt{ct>k,4>k) = 0. (95) 

Since the eigenvalues of pt are independent of t and since po is non-negative, 
Pt is non-negative as well. Hence \\pt\\i = Trpt and is independent of t. There- 
fore the local well-psedness of dH can be extended to the global one. 

The conservation of energy is proven in a standard way. □ 
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B Appendix: Jacobi identity for (071) 

Lemma B.l. The Poisson bracket defined in (I47p satisfies the Jacobi identity: 

{{A, B}, C} + {{C, A}, B} + {{B, C}, A} = 0. (96) 

Proof. In what follows we omit the argument p and the superindex c. Denote 
A' = dpA. Recall the definition, {A, B} = -jrTr {A' pB' - B'pA'). We have 

{A,B]='-Tr{[^,B']p). (97) 

We further denote A" = d'j^A. We think about A" as an operator on a two- 
particle space, or a two-particle operator. If we let Ki = K®I and K2 = I®K, 
then we have (omitting from now on the factor ^) 

{A^B}' = rr2(([A",i?^] + [4,5"])P2) + [A'^B'l (98) 

where Tr2 denotes the partial trace with respect to the second coordinate. The 
last two relations give 

{{A,B},C} = Tr{[{A,B]\C']p) 

= Tr,{[Tr2{[A", B^]p2 + [A, B"]p2), C[]pi) + Tr{[[A', B'], C']p). (99) 
Ler Tr*^^ denote the trace over the two-particle space. We have furthermore 

Tr,{[Tr2{[A",B'2]p2),C'M 

= Tr'^H[[A",B'2],C'Mp2) 
= Tr^\[[A",B[],C'2]piP2). 

and similarly for the second term on the r.h.s. of (j98p . Let R := p p and 
denote Tr®'^{Kp p) = TrR{K). Then we have, 

{{A, B], C} = TrniWr, B[IC'2] + [[A, B"l C[]) + Tr{[[^ , B'l C']p). 

Since the commutator of operators satisfies the Jacobi identity, the last equation 
implies that 

{{A, B], C} + {{C, A},B] + {{B, C}, A} 
= TrR{[[^\B[lC'2] + [[B",C[IA!2] + [[C\ J^^], B'2] 
+ [[A'^,B'%C[] + [[B'^^C'IA^] + [[C'2,A"IB[]) = 0. (100) 

□ 
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C Appendix. Proof of Eqn ( TMf ) 

Lemma C.l. 

Proof. We derive this inequality from the foUowing lemma: 
Lemma C.2. 

1 V (^-^)' V (^-^)' ^+-fiy no2) 

Proof. Let 

^-l.(Ar_p_n)!iV" UJ ^ ^ 

Then, 



(iV-p)! {N-p)l (1 

n=l ^ ' ^ ^ ?i=0 

(104) 

Subtracting equation (jlOSp from (1104p . we obtain that 

B = N-p_iN-pl Ay-^/y^"' (N-pV. N-p-n fl 

N iV^-P \2j ^ (A^-p-n)!iV"^ V2 

p (AT-p)! /n^-^' {N-p)l p + n fl 

~ ^ ~ iV ~ A^^-P V2 J ~ ^ (A^-p-n)!Ar"~/V" 

= l_Y^ {N-p)\ p + n(lY 



(N-p-n)\m N \2 

n=0 ^ ^ ' ^ 

which gives immediately equation (jl02p . □ 
On the other hand we have that 

N-p 



{N — p)\ p + nfl 



n=0 



{N-p-n)\m N \2 



N \2 - N 

n=0 



where in the last step we used equation (j42|) with p ^ p + \. This, together 
with equation (jl02p gives Lemma lC.ll □ 
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